Optimization by Box–Behnken design for environmental contaminants removal using magnetic nanocomposite

In this study, a CoO–Fe2O3/SiO2/TiO2 (CIST) nanocomposite was synthesized and utilized as an adsorbent to remove methylene blue (MB), malachite green (MG), and copper (Cu) from aqueous environments. The synthesized nanocomposite was characterized using field emission scanning electron microscopy (FE-SEM), Fourier-transform infrared spectroscopy (FTIR), thermogravimetric analysis (TGA), and X-ray diffraction (XRD). Input parameters included pH (3–10), contact time (10–30 min), adsorbent amount (0.01–0.03 g), and pollutant concentration (20–60 mg L−1). The effects of these parameters on the removal process efficiency were modeled and optimized using the response surface methodology (RSM) based on the Box–Behnken design (BBD). The RSM-BBD method demonstrated the capability to develop a second-degree polynomial model with high validity (R2 ˃ 0.99) for the removal process. The optimization results using the RSM-BBD method revealed a removal efficiency of 98.01%, 93.06%, and 88.26% for MB, MG, and Cu, respectively, under optimal conditions. These conditions were a pH of 6, contact time of 10 min, adsorbent amount of 0.025 g, and concentration of 20 mg L−1. The synthesized adsorbent was recovered through five consecutive adsorption–desorption cycles using hydrochloric acid. The results showed an approximately 12% reduction from the first to the seventh cycle. Also, MB, MG, and Cu removal from real water samples in optimal conditions was achieved in the range of 81.69–98.18%. This study demonstrates the potential use of CIST nanocomposite as an accessible and reusable option for removing MB, MG, and Cu pollutants from aquatic environments.


Materials and instrumentation
The materials used in the research including ethanol (C 2 H 6 O), cobalt nitrate hexahydrate (Co(NO 3 ) 2 •6H 2 O), hydrochloric acid (HCl), titanium (IV) butoxide (C 16 H 36 O 4 Ti), sodium hydroxide (NaOH), ammonia (NH 3 ), acetonitrile (C 2 H 3 N), methylene blue (C 16 H 18 ClN 3 S), iron (III) nitrate nonahydrate (Fe(NO 3 ) 2 •9H 2 O), potassium nitrate (KNO 3 ), tetraethyl orthosilicate (SiC 8 H 20 O 4 ), malachite green (C 23 H 25 ClN 2 ), and methanol (CH 4 O) were of laboratory purity and obtained from Merck, Germany.Stock solutions (1000 mg L −1 ) of each analyte in water were prepared through double distillation.These stock solutions were utilized to prepare the required concentrations for experiments and standard solutions.NaOH (0.1 M) and HCl (0.1 M) were used to adjust the pH of solutions, and pH measurement was conducted using a pH meter.For adequate mixing and contact between the adsorbent and analyte, a mixer with an intensity of 150 rpm at a temperature of 25 °C was employed.Centrifugation and a magnet were utilized to separate adsorbent particles from the solution during the reaction.The remaining concentration of dyes and metal ions in the samples was measured using UV-Vis spectrophotometry and atomic absorption spectroscopy, respectively.The characteristics of the adsorbent, such as morphology and phase features, were examined using FE-SEM, FTIR, TGA, and XRD.

Synthesis of CoO-Fe 2 O 3
Nano-absorbent was synthesized by hydrothermal method.First, 0.62 g of Co(NO 3 ) 2 •6H 2 O and 1.71 g of Fe(NO 3 ) 2 •9H 2 O were dissolved in 25 mL of double distilled water.Then the mixture was stirred on a stirrer for 1 h at room temperature.Next, NaOH 1 M was slowly added to the desired mixture to adjust its pH to 11.Again, the mixture was stirred for 1 h.Then the solution was transferred to the autoclave and placed under the temperature of 180 °C for 12 h.

Synthesis of CoO-Fe 2 O 3 /SiO 2
In order to synthesize CoO-Fe 2 O 3 /SiO 2 , 1 g of dried CoO-Fe 2 O 3 powder was mixed with 50 mL of deionized water, and then 150 mL of C 2 H 6 O and 2.5 mL of NH 3 were added to it.The mixture was dispersed under ultrasonic waves for 30 min.Then, while the solution was stirred under a mechanical stirrer, a mixture of 1.5 mL of SiC 8 H 20 O 4 and 0.5 mL of NH 3 was added dropwise to the solution.The reaction was continued for 6 h with continuous stirring at room temperature.Then, the formed particles were washed several times with C 2 H 6 O and

Experimental procedure
This study was conducted at the laboratory scale and in batch mode.The experiments were carried out based on the design provided by RSM.For this purpose, various concentrations of the analytes were prepared in several batches (250 mL).The pH of the samples was studied within the range of 4-10.HCl and NaOH were used to adjust the pH of the samples.In each adsorption experiment, a predetermined amount of adsorbent specified by the software was added to each batch and mixed on a shaker.After the required contact time, the adsorbent was separated from the solution by centrifugation (5 min, 2500 rpm) and using a magnet.Finally, samples were taken from each batch, and the concentration of dyes and metal ions was measured by UV-Vis spectrophotometry and atomic absorption spectroscopy, respectively.The percentage removal of analytes was calculated from Eq. 1.
where C 0 was the initial concentration (mg L −1 ), and C e was the concentration at the equilibrium time (mg L −1 ).

Characterization of the CIST nanocomposite
The morphology of CIST nanocomposite was examined through FE-SEM.Fig. S1a provides an image of the adsorbent surface, revealing nearly spherical and uniform magnetic nanocomposite with a particle size of approximately 20 nm, as evident in Fig. S1a.This uniform dispersion and the tiny size of nanoparticles contribute to the creation of a high surface area and a porous magnetic structure in the adsorbent.Additionally, the surface features, such as the present cavities, enhance its capability as an adsorbent.
Thermal stability is a crucial factor affecting the effectiveness of the adsorbent.Thermal analysis of CIST in the temperature range of 30-600 °C is presented in Fig. S1b.These nanocomposites exhibit thermal stability up to 600 °C.CIST shows a slight weight loss at 180 °C, attributed to the negligible weight reduction of metal oxides, with the minimal loss being associated with residual water molecules within the internal sheets and the surface of the nanomaterial.
The phase structure of CIST was examined through XRD analysis.The results are presented in Fig. S1c.In the diffraction pattern of CIST nanocomposite, distinctive peaks at 2θ values of 31.7°,36.2°,43.9°, 52.8°, 58.1°, and 65.2° indicate the crystalline structure of these magnetic nanocomposites.Peaks observed at 37.7°, 37.9°, 47.8°, and 54.3° correspond to the presence of titanium dioxide coating.Higher peaks at 2θ around 23.6° are also associated with the silica coating 31,32 .The particle size was determined to be 30 nm by Scherrer's equation.
The FT-IR spectrum of CIST in the wavelength range of 400-3400 cm −1 is illustrated in Fig. S1d.Spectrum analysis reveals a prominent peak at approximately 603.5 cm −1 , attributed to the presence of Co-O and Fe-O bonds.The peak observed at 1623 cm −1 is associated with water molecules in the inner layers of the nanoparticle structure, while the broad peak at 3284.6 cm −1 is indicative of stretching vibrations of O-H in silanol groups and water molecules.The absorption peak at 1099.5 cm −1 corresponds to the stretching vibration of Si-O-Si bonds, and peaks around 756.1 cm −1 indicate stretching and bending vibrations of Ti-O-Si, indicating the presence of titanium dioxide and silica coatings on the surface of CIST nanocomposite 33,34 .

Box-Behnken design
After selecting the effective factors, the BBD method was used to optimize the removal of MB, MG and Cu under the conditions of adsorbent amount, pH, analyte concentration and contact time.Each of the independent variables was defined at three levels: low (− 1), medium (0), and high (+ 1), based on coded values (Table 1).
The designed experiment included 27 different experiments, and the removal efficiency was considered as the response surface to the variables.In Table S1, the results of the conducted experiments using BBD, along www.nature.com/scientificreports/with the obtained removal percentages and predicted removal percentages, are presented.In the RSM approach, a model is defined for each dependent variable, expressing the main and interaction effects of factors on each variable separately.In this study, a second-order polynomial experimental model was used to investigate the four experimental variables.The response of the experimental system was based on Eq. 2.
In this equation, Y was removal efficiency, X i , and X j were the coded values of the factors, β 0 was the constant, and β i , β j , and β ij were the regression coefficients calculated by the software 35,36 .After selecting the appropriate mathematical model, statistical analysis of the data, examination of the adequacy of the experimental model, and optimization plot drawing were performed using the Design Expert software.Statistical data analysis was performed to assess the appropriateness of the applied experimental model.The ANOVA results for the polynomial response surface experimental model are presented in Tables 2, 3 and 4. As can be seen, the high quantities of the F-value and the obtained Probability ˃ F values (< 0.0001) denote the significance of the applied model 37 .Additionally, as shown in Tables 2, 3 and 4, the non-significant Lack of Fit mean squares indicate the adequacy of the second-degree model and the absence of other relationships affecting the removal efficiency.In conclusion, it can be inferred that the results obtained with the applied experimental model are in complete agreement.
As shown in Tables 2, 3 and 4, the regression coefficients for the linear fit R 2 were 99.87%, 99.43%, and 99.70% for MB, MG, and Cu, respectively.This indicates that over 99.4% of the variations in the removal are accounted for by the used dependent variables, and less than 2% of these variations cannot be justified by the experimental model.Moreover, the adjusted R 2 ˃ 0.98 for all three analytes was very close to R 2 , implying a good fit for the model.Therefore, the high values of R 2 and adjusted R 2 provide evidence for the appropriateness of the applied experimental model in this study.After ensuring the adequacy of the employed polynomial experimental model, multiple regression analysis using Design Expert software was conducted to determine the coefficients of the regression equation.Finally, solving the second-degree multivariable equation was employed to find the optimal (2)  In this equation, A was the amount of adsorbent (g), B was pH, C was concentration (mg L −1 ), and D was time (min).
In addition to the mentioned criteria for evaluating the accuracy of the selected model, the difference between predicted and experimental responses (residuals) has been graphically utilized to assess the model's precision.Residuals are considered as the deviations not captured by the model.The residual plots for evaluating the normal distribution of residuals are presented in Figs.S2a-c.The points in the residual normality plots form a straight line, confirming that residuals are normally distributed.Figs.S2d-f depict the predicted values against the actual values of the studied variables for MB, MG, and Cu.As observed, there is sufficient agreement between the actual data and the model-derived data, indicating the suitability of the selected model for the dataset.

3D response surface
The adsorbent mass is one of the influential factors in surface adsorption examined in the experiments.The results obtained from this study (Fig. 1a) indicate that by increasing the adsorbent mass from 0.01 to 0.03 g, the removal efficiency of MB increases.This enhancement in MB removal efficiency was observed at a contact time of 10 min and pH of 6 by varying the mass of CIST.The reason for the increased MB removal efficiency is that the higher adsorbent mass leads to an increase in the number of active adsorption sites in the solution.Consequently, the increased number of active sites enhances the contact surface between the adsorbent and the www.nature.com/scientificreports/pollutant, resulting in improved pollutant removal efficiency.Similar results were obtained in the studies by Qu et al. 38 and Kadhom et al. 39 .Another influential factor in the adsorption process is pH.This is because it induces a change in the ion state of the analyte, ionization, and surface charge of the adsorbent, affecting the interaction between the adsorbent and the adsorbate substance.In this study, it was observed (Fig. 1b) that the removal efficiency was low at very low pH levels.The reason for the low efficiency at acidic pH is creating a positive charge on the adsorbent, leading to an electrostatic repulsive force between the adsorbent and MB.Additionally, the concentration of hydrogen ions in the solution increases under acidic conditions, competing with MB for adsorption sites, resulting in decreased MB removal efficiency.However, as the pH increases, the concentration of hydrogen ions decreases, leading to increased MB adsorption.The study demonstrated that with an increase in pH, the removal efficiency of MB also increases.The removal efficiency of MB significantly increases by changing the pH from 4 to 6. Similar results were obtained in the studies by Jawad and Surip 40 and Mehmandost et al. 41 , indicating the similarity of the present study with their findings.
The required time for the interaction between the adsorbent and the adsorbate is a crucial factor that needs to be investigated in adsorption studies.The results of this study (Fig. 1c) demonstrated that with an increase in contact time, the removal efficiency of MG increases.This is attributed to the fact that, with an extended contact time, dye molecules have more opportunities to interact with the adsorbent surface.Initially, the adsorption capacity increases with time, but after a while, it reaches a relatively constant level, indicating the attainment of equilibrium.The initial increase in adsorption rate is due to the availability of more active sites and functional groups during the early stages.The study results revealed that the equilibrium time for this adsorbent is approximately 10 min.In a study conducted by Azam et al. in 42 on the removal of Cr and Cd using treated date seeds, it was found that the removal efficiency of Cr and Cd increases with contact time.Moreover, the highest adsorption occurs within the first 5 h.A similar outcome was observed in the study by Shojaei et al. in 43 , where they found an increase in removal efficiency with prolonged contact time until reaching equilibrium, which took approximately 20 min in their study.Afterward, no significant changes were observed.
Another factor investigated in this study was the pollutant concentration.The effect of Cu concentration on the removal efficiency was examined in the range of 20-60 mg L −1 , and the results are depicted in Fig. 1d.The study results revealed a decrease in removal efficiency with an increase in concentration from 20 to 60 mg L −1 .The reason for the reduction in removal efficiency at higher concentrations is that the active adsorption sites on the adsorbent's surface become saturated with the pollutant.As a result, the collision between analyte ions and active sites decreases.The findings of this part of the study align with the studies conducted by Mustafa et al. 44 and Gajera et al. 45 .

Optimization
In the removal process of MB, MG, and Cu, achieving the maximum removal was the primary objective in the experimental analyses.The optimal operational conditions were determined using numerical optimization techniques.Initially, optimization objectives, response surfaces, and independent variables were defined.The desirability function method was employed to obtain the best responses.The optimal removal conditions included an adsorbent dosage of 0.025 g, pH of 6, a concentration of 20 mg L −1 , and a contact time of 10 min.Under these conditions, the removal efficiency for MB, MG, and Cu was 98.01%, 93.06%, and 88.26%, respectively (Table 5).

Desorption studies
The appropriate eluent must be capable of thoroughly washing the adsorbed analyte from the adsorbent, allowing the adsorbent to be ready for reuse in experiments.In this study, tetrahydrofuran (THF), acetonitrile (ACN), methanol (MeOH), and hydrochloric acid (HCl) were employed as eluents.After each experiment, the adsorbent was separated from the solution using a magnet and placed in contact with the desired solvents to examine the desorption efficiency.The results of the eluent investigation on the analyte desorption efficiency are depicted in Fig. 2. According to the results in Fig. 2, HCl was suitable for desorbing MB, MG, and Cu from the surface of CIST, and it was utilized as the eluent in the experiments.

The determination of zeta potential
The point at which the positive and negative charges are equal is the zero point of charge (pH zpc ) for the adsorbent.At pH values higher than the pH zpc , the surface charge potential on the adsorbent is negative, and at pH values lower than the pH zpc , the surface charge potential is positive.To determine the pH zpc , 50 mL of 0.01 M KNO 3 was added to several beakers.Subsequently, the solution pH in the range of 2-10 was adjusted using 0.1 M HCl and 0.1 M NaOH.Next, 0.2 g of the adsorbent was added to each specified beaker and left to equilibrate for 24 h.The final solution pH was then measured using a pH meter, and a plot of initial pH versus final pH was generated.The intersection point of the two curves was identified as the pH zpc .The results of investigating www.nature.com/scientificreports/ the pH zpc for the CIST adsorbent are presented in Fig. S3.According to Fig. S3, the pH zpc for the adsorbent was determined to be 5.1.

The reusability of CIST nanocomposite
From a practical perspective, the recovery or disposal of pollutants from the adsorbent surface in the surface adsorption process is beneficial.The recovery and reuse of the adsorbent are crucial factors for evaluating its performance in practical applications.Therefore, after adsorbing MB, MG, and Cu from the aqueous solution with the CIST adsorbent, the adsorbent was recovered by washing with HCl and reused.The results of consecutive adsorption-desorption cycles for MB, MG, and Cu using the CIST adsorbent are shown in Fig. 3.According to the obtained results, the percentage of analyte desorption decreased by 12% until the seventh cycle.These results demonstrate the significant stability of the adsorbent and the non-destruction of the synthesized nanocomposite structure over several usage cycles.This indicates the economic viability of synthesizing and employing this adsorbent in treating wastewater containing MB, MG, and Cu.

Real samples analysis
To assess the efficiency of the proposed method in removing MB, MG, and Cu from real water samples, several water samples were utilized, including tap water, fish farm water, and wastewater.The procedure involved adding specified amounts of analytes to a fixed volume of water samples to prepare them for analysis.The desired samples were analyzed using the proposed method and under optimal conditions, with each experiment repeated five times.The obtained results are presented in Table 6.The results demonstrated that the removal efficiency of MB, MG, and Cu for real water samples ranged from 81.69 to 98.18%.The findings indicated that the sample matrix had a limited impact on the removal of MB, MG, and Cu from real water samples.

Comparison CIST nanocomposite with other adsorbents
According to the results of the present study, CIST nanocomposite can be used as an adsorbent for the effective removal of MB, MG, and Cu from water solutions.Therefore, the performance of CIST nanocomposite was compared with other adsorbents used to remove MB, MG, and Cu (Table 7).As shown in this table, CIST nanocomposite was better or comparable to other adsorbents in removing MB, MG, and Cu.Also, using RSM leads to a reduction in process time, a reduction in the number of tests, and a reduction in the use of materials.

Conclusion
In the present study, a CIST nanocomposite was synthesized and characterized by various techniques (XRD, TGA, FTIR, pH pzc , and SEM).Additionally, the influential factors affecting the removal, such as contact time, concentration, pH, and the amount of adsorbent, were investigated using RSM based on BBD.The removal efficiency of MB, MG, and Cu was examined under optimal conditions (i.e., pH of 6, contact time of 10 min, adsorbent amount of 0.025 g, and concentration of 20 mg L −1 ).The removal percentages of MB, MG, and Cu from real wastewater ranged from 81.69 to 98.18%.The reusability assessment of CIST nanocomposite demonstrated its suitable stability for repeated use in removing MB, MG, and Cu.According to the results obtained in this study, the synthesized nanocomposite can be considered a renewable adsorbent with significant impacts on wastewater treatment containing MB, MG, and Cu, making it an efficient material for removing these pollutants. https://doi.org/10.1038/s41598-024-57616-8

Figure 2 .Figure 3 .
Figure 2. The effect of eluent solvent on the desorption process.

Table 2 .
ANOVA for removal of MB.

Table 3 .
ANOVA for removal of MG.

Table 4 .
ANOVA for removal of Cu.

Table 6 .
Results of real samples under the optimum conditions (n = 5).

Table 7 .
Comparison of the CIST nanocomposite with other adsorbents.